Enhanced Screening in Chemically Functionalized Graphene 
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Resonant scatterers such as hydrogen adatoms can strongly enhance the low energy density of 
states in graphene. Here, we study the impact of these impurities on the electronic screening. We 
find a two-faced behavior: Kubo formula calculations reveal an increased dielectric function e upon 
creation of midgap states but no metallic divergence of the static e at small momentum transfer 
q — > 0. This bad metal behavior manifests also in the dynamic polarization function and can be 
directly measured by means of electron energy loss spectroscopy. A new length scale l c beyond 
which screening is suppressed emerges, which we identify with the Anderson localization length. 

PACS numbers: 72.80.Rj; 73.20.Hb; 73.61.Wp 
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Electronic screening presents a central problem in the 
physics of graphene and strongly affects electron trans- 
port as well as effects of electron-electron interactions in 
this material. First, external perturbations like charged 
impurities in the graphene substrate can be screened by 
the graphene electrons and the resulting scattering rate 
becomes inversely proportional to the squared effective 
dielectric constant e of graphene [H, 0. Second, electron- 
electron interactions are known to renormalize the charge 
carrier velocity in graphene near the Dirac point 0, Q 
and an excitonic instability can be realized depending on 

In reality, graphene samples are subject to disorder, 
which can alter e significantly. Single impurities can, for 
instance, lead to high energy plasmon poles of e @ . Most 
notably, impurities such as hydrogen adatoms can lead 
to pronounced resonances in the density of states near 
the Dirac point and related impurity bands can largely 
increase the low energy density of states (DOS) for impu- 
rity concentrations on the order of a few percent or less 
H&Eli- At high impurity coverages, however, graphene 
can turn into an insulating material such as graphane 
[ijj or fluorographene 13 with band gaps on the order 
of some eV 14| and thus vanishing DOS at low energies. 
Regarding the screening of electric fields in graphene it is 
unclear whether chemical functionalization with species 
like hydrogen metallizes graphene or rather turns it into 
an insulator. The answer to question is complicated by 
the special physics of Anderson localization in this ma- 
terial: The chiral symmetry has been shown to suppress 
Anderson localization in charge neutral graphene even in 
presence of strong local impurities (IR [lij . Away from 
the neutrality point, however, hydrogen adatoms can lead 
to Anderson localization 
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act Kubo formula calculations 0, [l9[. The dielectric 
function e(q,u>) and the dynamic polarization function 
Tl(q, ui) are discussed as function of momentum transfer 
q and frequency u>. As the central result, we find that res- 
onant impurities turn graphene into a bad metal: Screen- 
ing in graphene becomes indeed enhanced by chemi- 
cal functionalization, as the static dielectric function 
e(q,uj = 0) = e(q) exceeds the pristine graphene value 
of sc(q) ~ 4.93. Below a critical length scale l c , we find 
metallic behavior: e(q) <~ q^ 1 for q > q c = l/l c . How- 
ever, e(q) exhibits a maximum at q ~ q c and decreases 
with q for q < q c . Chemically functionalized graphene 
thus provides metallic screening only at length scales 
below l c . We further find bad metal characteristics in 
the dynamic polarization function and in — Iml/e(g,w) 
which can be measured by means of electron energy loss 
spectroscopy (EELS). By comparing to Green function 
calculations which neglect vertex corrections, we show 
that the bad metal behavior is not due to redistribution 
of single particle spectral weight but due to quantum in- 
terference processes. l c is identified with the Anderson 
localization length. 

To describe graphene functionalized with chemical 
species like hydrogen or organic groups, we consider the 
Hamiltonian H = H gr + Hi„ vp . H gr = — t J2<i j>( a l^j + 
h.c) is nearest-neighbor tight-binding Hamiltonian of 
graphene, where a\ (pi) creates (annihilates) an elec- 
tron on sublattice A (B) and t = 2.7eV is the near- 
est neighbor hopping parameter. The adsorbates are 
taken into account through the Hamiltonian Hi mp = 

e <iYl<iAdi + V J2i (d\ci + h.cj , where the parameters 
V = 2t, Cd = —t/16 for hydrogen have been deter- 
mined by density functional theory (DFT) calculations 
ll|. These impurities lead to resonances in the spec- 



In this letter, we study electronic screening in graphene 
with resonant scatterers by means of numerically ex- 



trum of graphene and to the formation of an impurity 
band, which is energetically very close to the Dirac point 
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Figure 1. (Color online) Static dielectric function of graphene 
with different concentrations of impurities m as function of 
wave vector q at fixed chemical potential p = 0. The maxima 
in e(q,uj — 0) indicate bad metal behavior. The wave vector 
q is in the Y — K direction. Periodic boundary conditions are 
used in a sample containing 4096 x 4096 carbon atoms. 



[llL IS, 2(| . Recently these impurity states have been 
observed experimentally by means of scanning tunneling 
spectroscopy (STM) [2l| and photoemission spectroscopy 
22|. 

The dynamical polarization function can be obtained 
from the Kubo formula [23| as 



n(q,w) 



dre^([p( q ,T),p(~ q ,0)}), (1) 



where A denotes the area of the unit cell, p (q) = 
J2i c l c i ex P ' r i) is the density operator, and the aver- 
age is taken over the canonical ensemble. For the case of 
a single-particle Hamiltonian Eq. ([I} can be calculated 
by means of the Chebyshev polynomial method 3, lif . 
The screening of electric fields is determined by the di- 
electric function £ (q, to). In the random phase approxi- 
mation (RPA) it reads e (q, ui) = 1 — V (q) U (q, lu) , where 
V (q) = 2ne 2 jq is the Fourier component of the Coulomb 
interaction in two dimensions. 

We start with analyzing the static dielectric function 
(Fig. [T]) obtained from our Kubo formula calculations. 
For pristine graphene, the RPA static dielectric function 
ec{q) = 4.93 is a constant [24[, which is also found in our 
calculations. In the presence of hydrogen impurities, one 
may expect that the screening is highly enhanced due to 
the emergence of midgap states in the vicinity of the neu- 
trality point and the correspondingly increased density of 
states around the Fermi level. Indeed, this expectation 
holds for impurity concentrations up to rii < 20% and 
the q-vectors (q > 0.05/a, where a is carbon-carbon dis- 
tance) accessible in our simulations 25|. However, the 



shape of the e((j)-curves does not show purely metallic 
behavior (i.e. e(g) ~ q^ 1 ), as might be expected from 
the non-zero density of states at the Fermi level: Wc 
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Figure 2. (Color online) (a,b) Static polarization function of 
graphene with different concentrations of impurities rii as a 
function of wave vector q. The chemical potential is fixed 
as n = 0. (c,d) g-dependent static polarization function of 
graphene at different chemical potentials p = 0, O.li, and 0.2t. 
Results for graphene without impurities (solid lines) are com- 
pared to those with an impurity concentration of rii = 0.5% 
(dashed lines) . (a,c) Results from the numerically exact Kubo 
formula calculations. (b,d) Results obtained by neglecting 
vertex corrections (Eq. ©)• 



find e(q) ~ q^ 1 but only for q > q c with impurity con- 
centration dependent q c . e(q) exhibits a maximum at 
q ~ q c and decreases for q < q c . Electronic screening in 
graphene is thus enhanced by hydrogen functionalization 
but fully metallic screening is found only at length scales 
below l c ~ l/q c . 

To understand the origin of this effect, we analyze the 
static polarization function (Il(q,u> = 0) = 11(g), Fig. [2]) 
obtained from the Kubo formula calculations and com- 
pare to the results obtained from evaluating a bubble 
diagram, 



n(q, iu n ) 



-y 

k,n 



TrG(k, icj n >)G(k + q,iu) n > +iw n ). 



(2) 

Here, G(k,iui n ) = (iuj n + fx — v F (a ■ k) - S(io; n )) 1 
is the Green function of graphene with Fermi velocity 
vp = 'Sat/ 2 in presence of strong local impurities, which 
are accounted for by the self-energy T,(i(j n ) = niT(iuj n ). 
In the limit of strong impurity potentials, the T-matrix 
entering the self-energy reads T(iui n ) = — I / G\ oc (iui n ) , 
where G® oc (iuj n ) is the local Green function of pristine 
graphene. This formalism accounts for spectral weight 
redistribution on a single particle level. The only approx- 
imation of Eq. ([2]) is that it neglects vertex corrections 
which describe quantum interference effects like Ander- 
son localization. 

While both approaches yield the similar H(q) for larger 
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Figure 3. (Color online) (a,b) Quasi eigenstates at energies 
E — — 0.031i (a) and +0.031i (b) for graphene with hydrogen 
adatoms at concentration n, = 0.5%. (c) Energy dependent 
localization lengths in graphene with different concentrations 
rii of hydrogen adatoms. (d) Critical screening length l c = 
l/q c , Anderson localization length £ and mean free path Le 
as function of average inter-impurity distance l/^/ru. l c is 
evaluated for /j, = 0, whereas £ and Le are evaluated for fj, 
inside the tails of the impurity band. 



q, we find that Eq. ([2]) yields constant non-zero H(q) for 
q — > for graphene with hydrogen impurities, which cor- 
responds to e (q = 0) — > oo. This is in contrast to finite 
and decreasing e (q) and decreasing Tl(q) for q — > as 
found by the Kubo formula simulations. The deviation of 
the bubble-diagram from the numerically exact Kubo for- 
mula results is strongest if the chemical potential (p « 0) 
is within the midgap impurity band generated by the hy- 
drogen adatoms (Fig. [2J; and d). Thus, the bad metal 
behavior of e(q) for q < q c must be due to a quantum in- 
terference effect involving the hydrogen induced midgap 
states. For doped graphene (/x ^ 0), the presence of the 
hydrogen impurities becomes less important. 

As Fig. [3] (a,b) illustrates, the eigenstates in the en- 
ergy range of the impurity band do not extend through 
the entire system but are localized in distinct regions of 
the sample. The spots where the quasi-cigenstates are 
localized depend strongly on energy, which is typical in 
the regime of Anderson localization [26[ . To corroborate 
Anderson localization as physical origin of the suppressed 
long wavelength screening near the neutrality point quan- 
titatively, we calculated electronic conductivities a, dif- 
fusion coefficients D(t, E) and the elastic mean free paths 
L e for undoped graphene (/i = 0) with different concen- 
trations of hydrogen impurities (see [III EH f° r techni- 
cal details). This allowed us to estimate the Anderson 



ization lengths (Fig. [3t): There is central peak showing 
comparably large localization lengths which paradoxially 
increase with impurity concentration. This peak is sur- 
rounded by plateaus, where the localization length is al- 
most independent of the energy. In the energy range of 
these plateaus, the localization lengths decrease with in- 
creasing impurity concentration as one expects. 

The comparison of the critical length scale for dielectric 
screening l c = l/q c to the Anderson localization length £ 
inside the plateau region (Fig. [3Ji) reveals a very similar 
dependence of l c and £ on the impurity concentration. 
It is thus very likely that the suppressed screening at 
long wavelengths is due to Anderson localization in the 
midgap impurity band. The critical length scale l c and 
the Anderson localization length correspond almost ex- 
actly to the average inter-impurity distance. It appears 
that this is the only natural length scale induced by the 
impurities. 

The increase of localization lengths with the impurity 
concentration in the center of the impurity band does 
not manifest in any of the dielectric screening properties 
investigated, here. It is plausible that the high DOS in 
the center of the impurity band leads to shorter average 
hopping distances between localized impurity states than 
in the tails of the impurity bands and might explain the 
increase of the conductivities / localization lengths with 
impurity concentration in that energy range. This mech- 
anism is qualitatively similar to variable range hopping 
}2(| with one difference. The latter requires some inelas- 
tic processes, usually electron-phonon scattering. Here, 
an energy uncertainty is provided by a general finite life- 
time broadening: The semiclassical conductivities a cor- 



respond by definition [151 . 1271 . 128j to the maximum con 
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localization lengths £ = L e exp(irh<j/2e 2 ) [if 
Inside the energy range of the midgap impurity states 
we find an interesting energy dependence of the local- 



ductivities obtained at finite simulation times and thus 
include an intrinsic lifetime broadening. 

The screening of static electric fields is determined by 
the generation of virtual electron hole pairs and, thus, 
involves processes from different energies. Being the real 
part of a retarded correlation function, Re II(q, ui = 0) 
can include information about virtual electron hole pairs 
at arbitrary excitation energies. Here, finite-energy par- 
ticle hole pairs contribute to Re n(q, uj = 0) due to fi- 
nite q, two C-p z bands being present in graphene and 
due to the disorder. Therefore, it is understandable why 
transport localization lengths evaluated at energies dif- 
ferent from the chemical potential in the screening calcu- 
lations display the same trend as l c in Fig. [3] Im n(q, ui) 
contains energy resolved information on processes con- 
tributing to the electrostatic screening and is closely re- 
lated to the electron energy loss function — Im (l/e(q, uj)), 
which can be measured by EELS. In normal metals 
we have ImII((j,g) ~ uj/q, while band insulators show 
Im H(w, q) = for energies ui less than the band gap. Un- 
doped graphene (fi = 0) lies between these cases: Here, 
only inter-band transitions are allowed, which leads to a 
peak at u> = v^q in the spectrum of Im n(q, a;). Below 
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Figure 4. (Color online) Dynamical polarization function 
of graphene at chemical potential (fj, = 0) with (a) differ- 
ent concentrations of impurities and the same wave vector 
q = 0.15a _1 as well as (b) with different wave vectors for the 
same concentration of impurities, (c) Electron energy loss 
function — Im(l/e) of graphene with adsorbed hydrogen im- 
purities at concentrations between m = and Hi = 2% with 
fixed q = 0.15a -1 . There is impurity induced EELS inten- 
sity at energies ui < VFq, which is energetically forbidden in 
pristine graphene. (d) — Im (1/e) in the low energy region for 
different momentum transfers q at fixed impurity concentra- 
tion 71; = 0.5%. 



the energy uj = vfq electron-hole excitations are forbid- 
den, as the energy exchange has to permit the momentum 
exchange hq. In the presence of hydrogen impurities, the 
midgap impurity band raises the possibility of electron- 
hole excitations in this forbidden region with spectral 
weight proportional to the impurity concentration m (see 
Fig. Ufa)). The peak at uj = v-pq is blue shifted and 
smeared out with larger impurity concentration. 

For fixed concentration of impurities, the slope of the 
low-frequency dependence of ImII(w,g) increases with 
decreasing q for q > q c (see Fig. 0Jb)) but remains con- 
stant with for q < q c . We have Imn(cj,g) ~ u/q for 
q > q c and ImII(u;, qj ~ u / q c for q < q c . Thus, the 
dynamical polarization function also shows bad metal 
characteristics with electronic-excitations being available 
at arbitrarily low energies but with non-metallic behav- 
ior at long wavelength/ small wave vectors q < q c . This 
bad metal behavior of chemically functionalizcd graphene 
is detectable in EELS experiments, as Fig. @] (c-d) 
show. The impurity band leads to an electron loss sig- 
nal —lm(l/e(q,ui)) > in the region ui < vpq which is 
energetically forbidden in pristine graphene. This loss 
signal increases with impurity concentration (c.f. 111(c)). 
As for ImLT(cj, q), we find metallic w-dependence, i.e. 
— Im(l/e(g,w)) ~ u for u ->■ (c.f. 01(d)). However, 
the momentum transfer dependence of the EELS signal 
is clearly non-metallic: we find — Im(l/e(g,a;)) increas- 



ing for q — > 0, whereas — Im (l/e(q,w)) is expected to 
approach a constant foro — >■ in the case of a normal 
two-dimensional metal [23 ] . 

In conclusion, we have studied electronic screening in 
graphene with resonant impurities. We show that the 
quasilocalizcd states in the vicinity of neutrality point 
lead to "bad metal" behavior. Resonant impurities make 
screening more efficient and render graphene metallic but 
only up to a certain critical length scale l c . This scale is 
determined by Anderson localization and corresponds to 
the inter-impurity spacing. As a consequence, scattering 
of the graphene electrons by external potentials, which 
vary on length scales which are large compared to the 
interatomic spacing but smaller than l c , is suppressed 
due to screening by the midgap states. 

Support by the Stichting Fundamcnteel Onderzock dcr 
Materie (FOM), the Netherlands National Computing 
Facilities foundation (NCF) and by DFG through SPP 
1459 are acknowledged. TOW thanks S. Kettcmann for 
useful discussions and RU Nijmegen for hospitality dur- 
ing the visit, when the ideas presented in this work were 
conceived. 



[1 
[2 

[3; 
[4 

[5 
[6 

[7: 

[8 
[9 

[10 
[11 

[12 
[13 



* |s.yu an@scien ce~ru.nl| 
1 |wehling@itp.uni-bremen.de| 
S. D. Sarma, S. Adam, E. H. Huang, and E. Rossi, Rev. 
Mod. Phys. 83, 407 (2011). 

M. I. Katsnelson, Graphene: Carbon in Two Dimensions 
(Cambridge Univ. Press, Cambridge, 2012). 
J. Gonzalez, F. Guinea, and M. Vozmediano, 
|Nucl. Phys. B 424, 595 (1994) | 

D. C. Elias, R. V. Gorbachev, A. S. Mayorov, S. V. Mo- 
rozov, A. A. Zhukov, P. Blake, L. A. Ponomarenko, I. V. 
Grigorieva, K. S. Novoselov, F. Guinea, and A. K. Geim, 
Nature Phys. 7, 701 (2011). 

J. E. Prut and T. A. Lahde, 
|Phys. Rev. Lett. 102, 026802 (20097] 
R. A. Muniz, H. P. Dahal, A. V. Balatsky, and S. Haas, 
|Phys. Rev. B 82, 081411 (2010) | 

V. M. Pereira, F. Guinea, TT M. B. Lopes dos San- 
tos, N. M. R. Peres, and A . H. Castro Neto, 
|Phys. Rev. Lett. 96, 036801 (2006)] 
N. M. R. Peres, F. Guinea, and A. H. Castro Neto, 



|Phys. Rev. B 73, 125411 (2006) | 

T. O. Wehling, A. V. Balatsky, M. I. Katsnelson, A. I. 
Lichtenstein, K. Scharnberg, and R. Wiesendanger, 
Phys. Rev. B 75, 125425 (2007). 

T. O. Wehling, M. I. Katsnelson, and A. I. Lichtenstein, 
Chem. Phys. Lett. 476, 125 (2009). 

T. O. Wehling, S. Yuan, A. I. Lichtenstein, A. K. Geim, 
and M. I. Katsnelson, Phys. Rev. Lett. 105, 056802 
(2010). 

D. C. Elias, R. R. Nair, T. M. G. Mohiuddin, S. V. 
Morozov, P. Blake, M. P. Halsall, A. C. Ferrari, D. W. 
Boukhvalov, M. I. Katsnelson, A. K. Geim, and K. S. 



Novoselov, Science 323, 610 (2009) 



R. R. Nair, W. Ren, R. Jalil, I. Riaz, V. G. Kravets, 



5 



L. Britnell, P. Blake, F. Schedin, A. S. Mayorov, 
S. Yuan, M. I. Katsnelson, H.-M. Cheng, W. Strupin- 
ski, L. G. Bulusheva, A. V. Okotrub, I. V. Grigorieva, 
A. N. Grigorenko, K. S. Novoselov, and A. K. Geim, 
|Small 6, 2877 (2010)1 

[14] M. Klintenberg, S. Lebegue, M. I. Katsnelson, and 
O. Eriksson, Phys. Rev. B 81, 085433 (2010). 

[15] F. Evers and A. D. Mirlin, 



|Rev. Mod. Phys. 80, 1355 (2008) 



[16] P. M. Ostrov sky, M. Titov, S. Bera, I. V. G ornyi, and 



A. D. Mirlin, |Phys. Rev. Lett. 105, 266803 ( 2010) 

[17] J. Bang and K. J. Chang, 

|Phys. Rev. B 81, 193412 (2010)| 
[18] S. Yuan, H. De Raedt, a nd M. I. Katsnelson, 

|Phys. Rev. B 82, 115448 (2010)| 
[19] S. Yuan, R. Roldan, and M. I. Katsnelson, Phys. Rev. 

B 84, 035439 (2011). 
[20] V. M. Pereir a, J. M. B. Lopes dos Santos, and A. H. 

Castro Neto, |Phys. Rev. B 77 115109 (2008) | 
[21] M. M. Ugeda, I. Brihuega, F. Guinea, and J. M. Gomez- 



Rodriguez, |Phys. Rev. Lett. 104, 096804 (2010) . 
[22] D. Habere^ L. Petaccia, M. Farjam, 37^ Taioli, 
S. A. Jafari, A. Nefedov, W. Zhang, L. Calliari, 



G. Scarduelli, B. Dora, D. V. Vyalikh, T. Pich- 
ler, C. Woll, D. Alfe, S. Simonucci, M. S. Dres- 
sclhaus, M. Knupfer, B. Buchner, and A. Griineis, 
|Phys. Rev. B 83, 165433 (2011) 



[23] R. Kubo, |J. Phys. Soc. Jp n. 12, 570 (1957) 
[24] " 

[25] 



B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New 
Journal of Physics 8, 318 (2006). 

Due to limitations of present computer power, it is nu- 
merically too expensive to get reliable values of the po- 
larization function for wave vectors q < 0.05/a. 
B. I. Shklovskii and A. L. Efros, Electronic Properties of 
Doped Semiconductors (Springer, Berlin, 1984). 



[26] 

[27] P. A. Lee and T. V. Ramakrishnan, 
|Rev. Mod. Phys. 57, 287 (1985)' 

[28] X Lherbier, ST M.-M. Dubois, X. Declerck, 
S. Roche, Y.-M. Niquet, an d J.-C. Charlier, 
|Phys. Rev. Lett. 106, 046803 (2011)1 

[29] A normal metal yields ReII(g — > 0,ui = 0) = k and 
ImII(g — s> 0, to — > 0) = £to/q with constants k and 
£. In RPA, we have e(q,uj) = 1 — (2ne 2 /q)Tl{q, u>) and 
consequently Kee(q,uj) 2> lme(q,uj) if nq 2> £<^>. At 
sufficiently small frequencies it is thus — Im l/e(q, lj) = 
Ime(g,ij)/(Ree(g,tj)) 2 ~ ujq° . 



V. V 1 



I 

<t 0.010 



0.008 



0.006 



Q 

CO 

I 

U 

o 

^ 0.002 
i 

0.000 
0.0 



0.004 



